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Abstract 

A new family of simple, analytic solutions of self-similar ly expanding fireballs is found for systems 
with ellipsoidal symmetry and a direction dependent, generalized Hubble flow. Gaussian, shell like 
or oscillating density profiles emerge for simple choices of an arbitrary scaling function. New, 
cylindrically or spherically symmetric as well as approximately one dimensional hydro dynamical 
solutions are obtained for various special choices of the initial conditions. 
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Introduction — Hydrodynamics is describing the local conservation of matter, momentum 
and energy. Due to this nature, hydrodynamical solutions are applied to a tremendous 
range of physical phenomena ranging from the stellar dynamics to the description of high 
energy collisions of heavy ions as well as collisions of elementary particles. Some of the most 
famous hydrodynamical solutions, like the Hubble flow of our Universe or the Bjorken flow in 
ultra-relativistic heavy ion physics have the properties of self-similarity and scale-invariance. 
Heavy ion collisions are known to create three dimensionally expanding systems. In case of 
non-central collisions, cylindrical symmetry is violated, but an ellipsoidal symmetry can be 
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well assumed to characterize the final state. The data motivated, spherical 
symmetric hydrodynamical parameterizations and/or solutions of refs. 

are generalized here to the case of such an ellipsoidal symmetry, providing new 
families of exact analytic hydrodynamical solutions. 

The new family of self-similar ellipsoidal solutions — The non-relativistic (NR) hydro- 
dynamical systems are specified by the continuity, Euler and energy equations: 

d t n + V(nv) = 0, (1) 
<9tV + (vV)v = -(Vp)/{mn), (2) 
<9 4 e + V(ev) = -pVv. (3) 

Here n denotes the particle number density, v stands for the NR flow velocity field, e for 
the energy density, p for the pressure and in the following the temperature field is denoted 
by T. These fields depend on the time t as well as on the coordinates r = (r x ,r y , r z ). We 
assume, for the sake of simplicity, the following equations of state, 

p = nT, e = up, (4) 

which close the set of equations for n, v and T. The NR ideal gas corresponds to k — §. 
The new family of exact analytic solutions of the hydrodynamical problem are given for 
arbitrary values of k > 0, m > 0. The hydro solutions are determined by the choice of a 
positive function of a non-negative, real variable T(s), corresponding to the (dimensionless) 
scaling function of the temperature. The scaling variable s is defined as 

s = J(2 + yi + ~Z2- ( 5 ) 
There the scale parameters depend on time, (X, Y, Z) = (X(t),Y(t), Z(t)). The ellipsoidal 
symmetry of the solutions is reflected by the ellipsoidal family of surfaces given by the 



s = so = const equation. The temperature and the density field depend on the coordinates 
(r x ,r y ,r z ) only through the scaling variable s. 

The new family of elliptically symmetric solutions of fireball hydrodynamics is given by 

V 

n(t,r) = n yu(s), (6) 

f s (X Y Z \ 
v (^ r ) = \j^r x ,-r y) -r z I , (7) 

T(t,r) = T Q (^) ' T(s), (8) 
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T(s) 1 \ 2T Jo T(u) r 
where the constant no is given by no = n(t Q , 0), the typical volume of the expanding system 
is V = XYZ, the initial volume being Vq = V(to), the dimensionless scaling function of the 
density profile is denoted by v(s), the constant To is defined by To = T(to, 0) and a constant 
of integration is denoted by Tj. The definitions of no and T correspond to the normalization 
u(s = 0) = 1 and T(s = 0) = 1. Initially, only one of the temperature and density profiles 
can be chosen as an arbitrary positive function, the equations of state relates the density 
and temperature profiles, resulting in the matching condition for the profile functions, as 
expressed by eq. ©• 

The equations of motion of the scale parameters are 

T / V \ 1/k 

XX = YY = ZZ = ^ (^2 J . (10) 

This time evolution of the radius parameters X, Y and Z is equivalent to the classical 
motion of a particle in a non-central potential, governed by the Hamiltonian 

where the canonical coordinates are (X, Y, Z) and the canonical momenta are (P x , P y ,P z ) = 
m(X,Y, Z). This generalizes earlier results 0, from k = 3/2 and T(s) = 1 to and to 
arbitrary T(s) > and k > 0. 

The conservation of energy by the classical Hamiltonian motion determines the physical 
meaning of the constant of integration T as the initial potential energy corresponding to 
the initial internal energy of the fireball: 

Etot = y(X 2 + Y 2 + Zl) + KT h (12) 



where the total energy is denoted by E tot , and the initial velocities are denoted by 
(X , Y , Z ). Due to the repulsive nature of the potential, the coordinates diverge for large 
values of t. The asymptotic velocities tend to constants 9( of (X as , ^as, Z as ), 

Etot = ^{Xl + Yl + Zl). (13) 

This completes the specification of the new family of solutions of fireball hydrodynamics 
with ellipsoidal symmetry. The form of the dimensionless scaling function T{s) can be 
chosen freely from among the positive functions of a non-negative variable. This freedom 
corresponds to a freedom in the specification of the initial conditions. Thus (uncountably) 
infinite new solutions of NR hydrodynamics are found. 

Self-similarity of the elliptic hydro of solutions — Each of these new hydrodynamical 
solutions is scale invariant: 

, , Xq Yq Z Q 
r = (r x — , Ty —,r z — ), (14) 

n(t,r)=n(t ,r')^§^), (15) 
v x (t,r) = ^ !•')-_,... (16) 

r(t,r) = r (( „r')(^) 1/ ". (17) 

Scale invariance of these solutions is equivalent to their self-similarity. The profile functions 
depend on time only through the scale parameters (X, Y, Z) and on the coordinates only 
through the scale parameter s. 

Limiting cases — The various physically interesting limiting cases fall into two classes. 
The first class of limiting cases corresponds to various additional symmetry properties im- 
posed on the scaling variable of eq. The second class of limiting cases corresponds to 
various choices of T(s), the scaling function of the temperature. The Lagrangian equations 
of motion for the scale parameters, eqs. (fTUj) as well as the shape of the flow velocity field, 
eq. (J7|), are the same for all the choices of T(s). Trivial prefactors are given in eqs. ()6I8|) . 
Hence these equations will not be repeated in the forthcoming discussion. We provide some 
physically interesting examples for the matching pair of density and temperature scaling 
functions (z/, T) that satisfy eq. (0). 

Spherically symmetric family of solutions — By assuming that initially all the scale 
parameters as well as their time derivatives are equal, X = Y = Z = R and X = Y = 



Z = R , the ellipsoidal family of solutions reduces to the spherical family of solutions of 
ref. [7] with a scale parameter X = Y = Z = R. The scaling variable and the equations of 
motion simplify to 

„2 rp / j-,3\ 1 / k 



h RR = Uw) ■ (18) 
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This family generalizes the Zimanyi-Bondorf-Garpman (ZGB) solution the spherical 
Gaussian solution jjj, and the Buda-Lund type of hydro solutions [3| to arbitrary scaling 
functions T{s) > and equation of state parameters k > 0. If i?o = 0, and the asymptotic 
velocity of the expansion is R 2 ^ = (u) 2 = Ti/m, and if k = 3/2, the equation of motion of 
the scale parameter is solved in a simple form 7} as R 2 = Rq + (u) 2 (t — t ) 2 . 

Cylindrically symmetric family of solutions — By imposing cylindrically symmetric initial 
conditions, X = Y = i? t0 , and X = Y = R t0 , one finds that X = Y = R t) as the equations 



of motion preserve cylindrical symmetry. Introducing r t = ^r 2 + r 2 , the scale parameter s 
and the equations of motion for the longitudinal and the transverse scales read as 

• = 3 + £ **= z *=£(iBf ■ (l9) 

These generalize the equations of motion for scales of the cylindrically symmetric, De - 
Garpman - Sperber - Bondorf - Zimanyi (DGSBZ) solution of ref. Q] to arbitrary < k ^ 3/2 
and to arbitrary T(s) > 0. 

One dimensional expansions — The equations of motion of parameters (X, Y, Z) has been 
studied for the case of k = 3/2 in ref. Although the expansion is generally 3 dimensional, 
a big initial compression in one of the directions (r 2 ) was shown to result in an effectively 
1 dimensional expansion in this direction, corresponding to Landau type initial conditions. 
In this case, an analytic solution for the variable Z is given in eqs. (23-25) of ref. [9(, and 
the conditions of validity of this approximation were given there by eqs. (26-28). These 
equations can be further simplified during the late stage of the expansion, when acceleration 
effects are small. In this limiting case, both eqs. (|lUj) and the conservation of energy are 
satisfied by the asymptotic solution: 

Y , l -mZ 2 a ~ ^, (20) 

(21) 
(22) 
(23) 



x a 


^ X , Y a ~ 


X(t) 


~ X + X a t, 


Y(t) 


^Y + Y a t, 


Z(t) 


~ z + z a t. 
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For simplicity, here we utilized to = . This asymptotic approximate solution is valid if the 



conditions of validity of the 1 dimensional expansion given in ref. 



9] 



are satisfied simulta- 



neously with the following constraints: Z <C Z a t, X ^> X t and Y ^> Y t. Alternatively 



m/(3Tj)Z <C t <C min(A /Ao, Y /Y ). Under these conditions all the initial internal en- 
ergy is converted into kinetic energy in the longitudinal direction, while the kinetic energy 
in the transverse components is conserved during the time evolution. For A = Y = 0, a 
one dimensional expansion is obtained. 

Various choices for the shape of the temperature scaling function T(s) generate interest- 
ing forms of the hydrodynamical solutions in all the ellipsoidal, cylindrical, spherical or 1 
dimensionally expanding classes. 

Gaussian solutions — The simplest possible choice for the temperature scaling function 
is T(s) = 1. After a trivial scale transformation, (A, Y, Z) — > y^(A, Y, Z), 

u(s) = exp(-s/2), T(s) = l. (24) 

The density profiles are Gaussians and the temperature distribution becomes spatially homo- 
geneous, as follows from eq. ©, and we recover the elliptic Gaussian solutions described in 

nn 

refs. |9|,llQl|- If freeze-out happens at a constant value of the local temperature, T(t, r) = Tf, 
these Gaussian hydrodynamical solutions corresponds to a sudden freeze-out at a constant 
time, t — tf, in the whole volume of the fireball. Remarkable features of this model are that 
i) the slope parameters of the transverse momentum spectra increase linearly with mass, 
with the coefficient of linearity depending on the relative direction to that of the impact pa- 
rameter, ii) the parameters of the two-particle Bose-Einstein correlation functions oscillate 
as a function of the angle between the event plane and the transverse momentum of the 
pair |l0| . 

De - Garpman - Sperber - Bondorf - Zimdnyi solution — One can require that the 
temperature and the density profiles are described as different powers of the same profile 
function. Such a similarity is achieved if the temperature and the corresponding density 
profiles are 

T{s) = (1- a) 8(1 -a), (25) 
u(s) = (1 - s) a 8(1 - a), a = T,/(2T ) - 1. (26) 

These profile functions generalize the spherical ZGB solution Q], and the cylindrically 
symmetric DGSBZ solution to asymmetric ellipsoids and arbitrary values of k. The 



Elliptic fireball 




FIG. 1: An approximately Gaussian BL-H fireball profile is shown for X = 5 fm, Y = 8 fm, 
m/To = 1, (u) = 0.7, (AT/T) = 0.1. The vertical scale is arbitrary, the longitudinal coordinate is 
r z = 0. For the time dependence see eqs. (TJOJ) . 

physical meaning of the parameter a is determined here, in terms of the total initial internal 
energy Tj and the initial central value To of the temperature field. 

Elliptic Buda - Lund solutions — The Buda - Lund hydro model (BL-H) was developed 
for the description of single-particle spectra and two-particle Bose-Einstein correlation func- 
tions in high energy heavy ion collisions at CERN SPS 0, |^. The BL-H attempted to 
characterize the temperature, density and flow fields by their means and variances only, its 
essential property is that it simultaneously involves a temperature gradient parameter and 
a Hubble-like flow profile. In its original form, the model was cylindrically symmetric and 
the (longitudinal) flow profile was relativistic BL-H type of exact hydro solutions 

correspond to the scaling functions 

1 



T(s) 



1 + 

(1 + bs) exp 



1 AT 

2 \ T h 



T; 



2T 



[s + bs 2 /2) 



(27) 
(28) 



The dimensionless parameter b is interpreted as a measure of the transverse temperature 
inhomogeneity Q, 11|. Regardless of the symmetry classes, the BL-H density profile has an 
approximately Gaussian, conventional shape if b < Tj/(2To). On the other hand, the density 
profile looks like an ellipsoidal ring of fire, with a density minimum at the center and a density 
pile- up on the surface, if b > T^/Tq. In the spherically symmetric case, similar morphological 
classes of BL-H solutions were found in ref Introducing the notation Tj = m(u) 2 , one finds 



that 



m{«} 2 



> (—) yields ellipsoidal, expanding fireball BL-H solutions, Fig. [T] A detailed 
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Elliptic shell of fire 




FIG. 2: A BL-H shell profile with X = 5 fm, Y = 8 fm, m/T = 1, (u) = 0.5, (AT/T) = 0.71, 
otherwise as Fig. 1. 

analysis of correlations and spectra in Pb + Pb collisions at CERN SPS indicated such 
a behavior, corresponding to big and expanding fireballs. On the other hand, the analysis 
of correlations and spectra of h + p reactions at CERN SPS indicated 12, ^3, Q] small 



transverse flow and big transverse temperature inhomogeneity, m ^ < (^r). This case 
corresponds to the formation of a shell of fire. Such an example is shown in Fig. |21 

Fireballs with density waves — Here we show that a simple choice of the temperature 
scaling function can lead to periodically oscillating temperature and density waves in exact 
solutions of NR hydrodynamics. A pair of oscillating scaling functions is, for example 

(1 + a cos Ps)^ 1 , 
(1 + a cos f3s) exp 



T(s) 



(29) 



V[8) 



-—^r(s — t: sin 3s) 
27V P 



(30) 



where parameters (a,/?) correspond to the amplitude and the period of the oscillations, 
respectively, as shown in Fig. 3. As time evolves, the shapes of the ellipsoids change 
in coordinate space: the bigger the initial compression, the faster the expansion in that 
direction, corresponding to eqs. (jlOj) . 

Connection with other exact solutions of fireball hydrodynamics - Before summarizing 
the results let us also relate the presented family of exact solutions of hydrodynamics to 
presently known, exact solutions of relativistic hydrodynamics. One of the most well known 
relativistic hydrodynamical solution is the 1+1 dimensional Landau-Khalatnikov solution, 
described in refs. 119] , Q and 0- 

The initial condition here is a uniformly heated piece 




FIG. 3: An oscillating wave-like density profile is obtained, for X = 5 fm, Y = 8 fm, Tj/To = 1, 
a = 0.2, = 2.0. 

of matter initially at rest, the equation of state is that of an ultrarelativistic gas with three 
degrees of freedom, but expanding only 1 + 1 dimensions. Landau has shown that such ex- 
pansions lead to an approximately Gaussian rapidity distribution. This Landau-Khalatnikov 
solution is presently the only known exact solution of relativistic hydrodynamics that de- 
scribes an exploding fireball with relativistic acceleration. These solutions are, however, not 
self-similar, and not explicit, one dimensional solutions, hence they are not easily related 
to the exact, explicit, selfsimilar and accelerating non- relativistic solutions described in the 
present manuscript. 

However, it is interesting to note, that in the late time limit, the acceleration of the scale 
parameters vanishes in the non-relativisitic self-similar solutions described here, as given by 
eqs. f)21H23j) . Hence for very late time, the velocity field becomes spherically symmetric, 
v — > r/t, which corresponds, within the lightcone, to a well known solution of relativistic 
hydrodynamics: u^(x) = x^/r, which is in the 1+1 dimensional case the reknowned vel ocity 
distribution of the Hwa-Bjorken solution of relativistic hydrodynamics, given in refs. [tsI]. |l q(| . 
while in 1 + 3 dimensions v = r/t or u^(x) = x^/t is referred to as the Hubble flow, as it 
corresponds to the flow velocity field of galaxies in an expanding Friedmann universe. These 
solutions are boost-invariant and both in the 1+1 dimensional Hwa-Bjorken and in the 1+3 
dimensional Hubble case, corresponding to a flat rapidity distribution. The initial boundary 
conditions are given only within a lightcone, on a r = r boost-invariant hypersurface, the 
equation of state is characterized by a constant or piecewise constant speed of sound, and 
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the solutions are explicit, accelerationless and self-similar flows. 

In ref. [2Q] these boost-invariant self-similar solutions were generalized to 1+1 dimen- 
sional, non-boost invariant, self-similar, accelerationless solutions, using a broad class of 
equations of state: e = mn + Kp, p = nT . The initial conditions were given, similarty 
to the case of the Hwa-Bjorken solution, on a boost-invariant hypersurface with r = r , 
but assuming an inhomogeneous initial temperature profile and a corresponding matching 
initial density profile. In the same paper this solution is extended to axially symmetic, 
three dimensionally expanding fireballs, that represent the late stages of central heavy ion 
collisions. However, for non-central collisions, axial symmetry is too restricitve and due to 
this reason we have generalized these solutions for ellipsoidally expanding relativistic fire- 
balls, and the 1+1 dimensional solution, the 1+3 dimensional axially symmetric solution 
and the 1+3 dimensional ellipsoidal solutions were written up also in a series of papers 
for the 3rd Budapest Winter School on Heavy Ion Physics, published in refs. |21|], 22], 23J. 
These self-similar, relativistic hydrodynamical solutions, in particular, the ellipsoidal expan- 
sion described in ref. 2^, correspond exactly to the late time limit of the non- relativistic 
solutions described here. See ref. J3] for further details on such a correspondence as well as 
for a more complete list of citations on early papers on the Hwa-Bjorken solution. 

It is worthwhile to mention that there are additional, recently found exact solutions of 
relativistic hydrodynamics that find accelerationless generalized Hubble type of solutions, 
with direction dependent Hubble constants. The first of these class of solutions has been 
found by Biro in refs. [2^, 2sl |. for cylindrically symmetric expansions at the softest point 
of the equation of state. Sinyukov and Karpenko recently published a solution which can 



be considered as the generalization of refs. 



2j,0, as well as that of ref. H, for the 



pre-asymptotic, accelerationless stage of the expansion, when the scales expand already 
linearly in time, corresponding to eqs. f)21H23j) . but before the time period when the off-sets 
(X a , Y a , Z a ) can be considered negligibly small. The Biro as well as the Sinyukov-Karpenko 
solutions are also self-similar and accelerationless relativistic solutions, their shortcoming 
is that they are valid in the medium only when the pressure is a temperature independent 
constant, or, when the boundary condition is an expansion to the vacuum, they are obtained 
only for the special "dust" equation of state, that corresponds to a uniformly vanishing 
pressure, p = p = 0. These solutions can thus be considered as a pre-asymptotic, relativistic 
solutions corresponding to very late stages of the fireball expansions. 
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This brief review of known exact solutions of fireball hydrodynamics in the relativistic 
kinematic domain indicates, that we find a correspondence between the late stages of the 
presented exact, accelerating, self-similar solutions of non-relativistic, ellipsoidally symmet- 
ric fireballs and the known class of relativistic but accelerationless, self-similar, ellipsoidally 
symmetric solutions of hydrodynamics. The missing link between them is the class of rel- 
ativistic, accelerating, explicit solutions of hydrodynamics that use a realistic equation of 
state. The search for this missing class of exact and explicit solutions has been started and 
the first results will soon be reported elsewhere. 

Summary — The non - relativistic hydrodynamical problem has been solved for expanding 
fireballs with ellipsoidal symmetry for the class of self-similar expansions. The flow velocity 
distribution is a generalized Hubble field in all the cases. An exact solution is assigned 
to each positive, integrable function of a non-negative variable. The time evolution of the 
(X, Y, Z) scale parameters corresponds to a Hamiltonian motion of a mass point in a non- 
central, repulsive potential. The density profiles may be of fireball type, or they may form 
one or more shells of fire. For initial conditions with higher symmetry, one dimensional, 
cylindrical and spherical expansions are obtained. 

These results provide analytic insight into the time evolution of expanding fireballs with 
nontrivial, ellipsoidally symmetric morphology, a non-polynomially hard algorithmic prob- 
lem that is difficult to solve even numerically on von Neumann type computers. 

Acknowledgments: Thanks are due to L. P. Csernai, Y. Hama, F. Grassi, M. Gyulassy, 
T. Kodama, B. Lukacs and J. Zimanyi for inspiring discussions. This work was supported 
by OTKA grants T026435, T034296, T 038406, NWO - OTKA grant N 25487, by an NSF - 
MTA - OTKA grant INT0089462, the US DOE grant DE - FG02 - 93ER40764 and by the 
NATO PST.CLG.980086 grant. 



* Electronic address: csorgo@sunserv.kfki.hu 
[1] J. Bondorf, S. Garpman and J. Zimanyi, Nucl. Phys. A296 (1978) 320 . 
[2] J.N. De, S.I.A. Garpman, D. Sperber, J.P. Bondorf and J. Zimanyi, Nucl. Phys. A305 (1978) 
226. 

[3] T. Csorgo, B. Lorstad and J. Zimanyi; Phys. Lett. B338 (1994) 134. 



11 



T. Csorgo and B. Lorstad, Phys. Rev. C54 (1996) 1390. 
T. Csorgo and B. Lorstad, Nucl. Phys. A590 (1995) 465c. 

J. Helgesson, T. Csorgo, M. Asakawa and B. Lorstad, Phys. Rev. C56 (1997) 2626. 

T. Csorgo, Centr. Eur. J. Phys. 2 (2004) 556-565 

P. Csizmadia, T. Csorgo and B. Lukacs, Phys. Lett. B 443 (1998) 21. 

S.V. Akkelin, T. Csorgo, B. Lukacs, Yu. M. Sinyukov and M. Weiner, Phys. Lett. B 505 
(2001), 64. 

T. Csorgo, S.V. Akkelin, Y. Hama, B. Lukacs, Yu. M. Sinyukov, Phys. Rev. C 67 (2003) 
034904 

A. Ster and T. Csorgo, Nucl.Phys. A661 (1999) 419-422 
T. Csorgo, |heFph/00012331 Heavy Ion Phys. 15(2002)1 



W. Kittel, Acta Phys. Polon. B32 (2001) 3927-3972 

N. M. Agababyan, EHS/NA22 Collaboration, Phys. Lett. B422 (1998) 359 

L. D. Landau, Izv. Akad. Nauk Ser. Fiz. 17 (1953) 51. 

I.M. Khalatnikov, Zhur. Eksp. Teor. Fiz. 27 (1954) 529; see also Q. 

S. Z. Belenkij and L. D. Landau, Nuovo Cim. Suppl. 3S10 (1956) 15 [Usp. Fiz. Nauk 56 
(1955) 309]. 

R. C. Hwa, Phys. Rev. D 10 (1974) 2260. 
J. D. Bjorken, Phys. Rev. D 27 (1983) 140. 

T. Csorgo, F. Grassi, Y. Hama and T. Kodama, Phys. Lett. B 565 (2003) 107 
|arXiv:nucl-th/0305059| . 

T. Csorgo, F. Grassi, Y. Hama and T. Kodama, Heavy Ion Phys. A 21 (2004) 53 [Acta Phys. 
Hung. A 21 (2004) 53] |arXiv:hep-ph/0203204| . 

T. Csorgo, F. Grassi, Y. Hama and T. Kodama, Heavy Ion Phys. A 21 (2004) 63 [Acta Phys. 
Hung. A 21 (2004) 63] |arXiv:hep-ph/0204300| . 

T. Csorgo, L. P. Csernai, Y. Hama and T. Kodama, Heavy Ion Phys. A 21 (2004) 73 
|arXiv:nucl-th/0306004| . 

T. S. Biro, Phys. Lett. B 474 (2000) 21 j|arXiv:nucl-th/9911004| . 



T. S. Biro, Phys. Lett. B 487 (2000) 133 |arXiv:nucl- th/0003027 . 



Y. M. Sinyukov and I. A. Karpenko, arXiv:nucl-th/0506002 



12 



